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It may well happen that the two main components of the dark sector of the Universe, dark matter and 
dark energy, do not evolve separately but interact nongravitationally with one another. However, 
given our current lack of knowledge on the microscopic nature of these two components there is 
no clear theoretical path to determine their interaction. Yet, over the years, phenomenological 
interaction terms have been proposed on mathematical simplicity and heuristic arguments. In 
this paper, based on the likely evolution of the ratio between the energy densities of these dark 
components, we lay down reasonable criteria to obtain phenomenological, useful, expressions of the 
said term independent of any gravity theory. We illustrate this with different proposals which seem 
compatible with the known evolution of the Universe at the background level. Likewise, we show 
that two possible degeneracies with noninteracting models are only apparent as they can be readily 
broken at the background level. Further, we analyze some interaction terms that appear in the 
literature. 

PACS numbers: 


PACS: 98.80.-k, 95.35.+d, 95.36.+X 


I. INTRODUCTION 

Usually the present Universe is described by the Friedmann-Robertson-Walker (FRW) metric coupled to the 
dominant energy components, viz. baryons, dark matter (DM), and dark energy (DE). The two first constitute the 
main ingredients of cosmic structure. The latter two comprise the so-called “dark sector” as their presence is felt 
only by their gravitational interaction with the former. While DM (alongside the baryons) is nearly pressureless, 
DE is endowed with a high negative pressure responsible for the present accelerated expansion. Since, thus far, 
the information available on these two components is indirect —they have never been directly detected, much less 
produced in the lab — their nature remains unveiled. This may explain the very often minimalist assumption 
that they evolve independently of each other, i.e., that they interact with one another (and with whatever other 
component) gravitationally only. Nevertheless, there is no reason a priori why they should absolutely not interact 
otherwise; further, a coupling in the dark sector is not only natural, it seems inevitable Q-i- 

The simplest cosmological model with a high degree of success from the observational side [the so-called lambda cold 
dark matter (ACDM)] takes as DE component the cosmological constant, or what is dynamically equivalent, the 
energy density of the quantum vacuum though finely tuned to a tiny value. Alongside this tuning problem it also 
presents the so-called “coincidence problem”, i.e., the fact that nowadays the energy density of DM is of the same 
order as that of the cosmological constant despite the former declines with cosmic expansion as a~^ (where a is the 
scale factor of the FRW metric) while the latter remains fixed Q . 

A way to alleviate, or even solve, the latter problem is to allow for some nongravitational interaction in the dark 
sector that gives rise to a continuous transfer of energy from DE to DM. (Note that a transfer in the opposite 
direction would only exacerbate the problem, so it is more appealing to assume that the transfer occurs just in 
that sense; although the opposite possibility is not excluded, neither is the possibility of one or more changes in the 
direction of the transfer as the Universe expands). Current data are compatible with such transfer although the 
evidence is, so far, not fully conclusive Eor reviews on the subject see [T^, [ll| and references therein. See, also. 
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the recent contributions 

Given our ignorance about the nature of the components of the dark sector (to begin with, we do not even know 
whether these are single or multiple components), there is no reliable guidance to determine the interaction term. 
Therefore, the expressions proposed in the literature (except for the obvious restriction that such term must be small, 
not to deviate seriously from the predictions of the ACDM model) are necessarily phenomenological and based on 
heuristic arguments. Here we propose a different route to arrive at reasonable expressions of such term based on the 
likely evolution of the ratio between the densities of DM and DE. 

From the history of structure formation we know that the former dominated at early times (though preceded by a 
short period of radiation prevalence), and from the current accelerated expansion we infer that the latter started to 
dominate in the recent past, and it will likely dominated for ever —at least if the ACDM or a slight variant of it turns 
to be the right cosmological model. Therefore, qualitatively speaking, the said ratio must decrease monotonously for 
most of the time. Although, as argued below, at very high redshifts it must have varied little or even stayed constant, 
and it will also scarcely vary or be constant in the far future. Accordingly, by wisely parametrizing the said ratio and 
using the conservation equations of the energy components, valid in most theories of gravity, compelling expressions 
for the interaction term are readily obtained. 

This paper is organized as follows. Section II relates the interaction term with the said ratio, establishes the criteria 
that any reasonable parametrization of the latter should comply with and illustrates it by proposing three different 
parametrizations based solely on the FRW metric and the energy conservation equations. Section III makes use of 
recent data on the Hubble history to constrain the free parameters entering these parametrizations. Section IV deals 
with two apparent degeneracies between the interaction term and the equation of state (EOS) of DE. Section V 
considers different proposals for the interaction term made in the literature to see whether their corresponding ratios 
between the energy densities comply with the criteria of SEc. H. Finally, Sec. VI summarizes our findings. 

As usual, a subscript zero attached to any quantity means that the latter is to be evaluated at present time. 


II. PARAMETRIZING THE RATIO DM/DE 

Let us consider a homogeneous and isotropic universe whose main energy components are baryons, dark matter and 
dark energy, subscripts b, m, and x, respectively. The first two being pressureless, the latter with a not necessarily 
constant EOS parameter w =Px/px < —1/3. 

We first look for an equation describing the dependence of the quantity, Q, that gauges the strength of the interaction 
in the dark sector, on the ratio between the densities of DM and DE, r = Pm!Px- Thus, we write the conservation 
equations of the components as 


Pb + 

3Hpb = 0, 

(1) 

Pm T 

pjn — Q ^ 

(2) 

Px + 

3H(1 -1- w)px = -Q , 

(3) 


where H = a/a denotes the Hubble function. In writing last two equations we have implicitly assumed that all DM 
and all DE participate in the interaction. For Q > 0 energy flows from DE to DM and vice versa for Q featuring the 
opposite sign; however, as mentioned above, we will take Q > 0 as it alleviates the coincidence problem. In accor¬ 
dance with strong constraints from local gravity experiments [a in we assumed that the baryons conserve separately. 


As it is readily derived, the condition Q > 0 imposes f/r > 3Hw, i.e., it restricts the rate at which r can decrease as 
the Universe expands. (Recall that in the ACDM model {w = —I) one has r = —3Hr). 


From the above, we get 


Q 



H 


Pm 

1 + r ’ 


(4) 


where the prime stands for d/da. 
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We reasonably expect Q to be smaller than the second term on the left-hand side of (2). This, alongside the 
restriction Q > 0, leads to 0 < ar' — Swr < 3(r -|- r^). 

We now look for reasonable expressions for r to be inserted in this, by passing, will illustrate our method. 
Any expression for r should satisfy dr/da < 0, except when a —> 0 and when a —> oo, because in these limits r 
must either state constant or vary very slowly. Hence the quantities r+ = r(a —>■ 0), ro = r{a = 1) ~ 25/70 and 
r_ = r(a —>• oo) must fulfill 0 < r_ < ro < 1 < r+, with r+ finite. Should the latter diverge, it would be because 
of one of the three following possibilities: (i) px vanished and pm stayed finite as a ^ 0. Then, the initial amount of 
DE would be zero and no transfer of DE to DM could ever occur whence Q would vanish altogether at all times, 
(ii) Px remained bounded but Pm diverged (in the same limit as above). In this second case Eq. (2) implies that Q 
would diverge, and by (3) Px would become minus infinity. This, in turn, would lead to the instantaneous transfer 
of the total amount of DE to DM. That is to say, the dark energy component would disappear completely for good. 
However, if pm would obey pm = Constanta”^ -I- e(a) where e does not vanish as a —>■ 0, things would be different, 
because in this instance, Q does not diverge. This introduces an extra uncoupled DM parameter that could be 
included when comparing with the data. Nevertheless, the expression for pm corresponds to the case where not all 
DM interacts, only a part of the total, but this was already excluded from our assumptions above, (iii) Both energy 
densities diverge, but pm tends to a higher infinity than px does. Clearly, this case is the same as the previous one. 
Einally, the expression for r must gentl y in terpolate between and r_. Observational data suggest that the former 
is much larger than unity; see, e.g., [la, [l7l|. 

The possibility that for some period between and rg the ratio r is not decreasing would imply a severe deviation from 
the ACDM model that, most likely, would bring observable consequences on structure formation and the integrated 
Sachs-Wolfe effect. This is why we assume that dr/da < 0, except possibly at r+ and r_. 


A. First parametrization 


A possible choice for the ratio Pm!Px is 


r = r. 


r+ — r_ 

1 -h /3a“ ’ 


(5) 


where a and /3 are positive-definite constants, 
and the corresponding current ratio is rg = 
w = constant), r+ ^ r_, then r+ ~ rg (1 -|- /3). 


As it can be readily checked ([5]) complies with the above conditions, 
r_ -I- (r+ — r_)(l -|- /3)“^. If, as it seems reasonable (at least for 
This suggests (3^1. 


From the condition Q > 0, it follows that ar'/r > 3w. This alongside the observationally well-supported assumption 
r+ ^ r_, yields 


r 


a/3a“ 

1 -h /3a“ 


> 3w , 


thereby in the limit a 1, the constraint a < 3 | w | follows. 

The first derivative of r, evaluated at present time (a = 1) and under the assumption r_ <C r+, yields 


( 6 ) 


r'(a 



a ~ —rga. 


(7) 


Because 0 < a < 3, it is seen that | r'(a = 1) |< 3rg whereby the slope of r at a = 1 is more gentle than the 
corresponding one of ACDM model, 3ro, whence the coincidence problem gets alleviated. 

Eigure[l] shows the evolution of r for some chosen values of the free parameters compatible with observational data, 
at least, at the background level. 

In the particular but interesting case that ic = — 1 and r+ r_ we have 


Q 



al3a°‘ \ 
1 + Pa°‘) 


l-F/3a“ 
r+ -|- 1 -|- /3 a“ 


( 8 ) 


^ For simplicity, we take the current energy budget as: baryons 5%, DM 25%, and DE 70%. 
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Figure 1: The first curve (from left to right) shows the evolution of ratio pm!Px corresponding to the expression ([5} where 
we have used a = 2.70 and p = 89.6, consequently r+ = 32.36.The same can be said for the next curve but with a = 2.76, 
P = 279, and r+ = 100. The right-most curve depicts for comparison, the evolution of the ratio for the concordance ACDM 
model, pm/pA = roa~^. In all three cases, we have taken ro = 25/70. As it is apparent, the interaction Q > 0 mitigates the 
coincidence problem. 


This illustrates how a reasonable expression for the interaction term between DM and DE can be devised from simple 
arguments based on the evolution of the ratio Pm!Px- 

By inserting the last expression in Eq. (2) we obtain, after integration, the evolution of the DM density, 

Q(l + r_|_) —3r_|_ 

Pm -4^ \l + r++Pa°‘ 

PmO L 1 + r_|_ -|- p 

In the limiting cases a <C 1 and a oo, one has pm oc and pm oc , respectively. These results 

are very reasonable. At early times the interaction has not had much chance to be felt; besides, the second term 
on the left-hand side of (2) is much larger than Q, thereby the DM density goes down practically as if it were no 
interaction, pm a~^. At late times, the interaction has had plenty of time to be felt and Q is not so small as 
compared with the said second term; thus pm departs from the case of just pure expansion and goes down more slowly. 

Figure [5] depicts, for illustrative purposes, the evolution of the ratios Pm/Pmo and Px!Pmo for some selected values of 
the free parameters. 


»(i+r+) 


(9) 


B. Second parametrization 


Another possible expression for the ratio Pm!Px is 


r = T- exp 



( 10 ) 


where r_ = r{a —>■ oo) > 0 and A > 0. On the other hand, r+ = r(a —>■ 0) = rp see Fig. [3] for the cases A = 8 
and 10, both with ro = 25/70. 

The condition for alleviating the coincidence problem | r'{a = 1) |< 3ro implies A < 12. For w = —1, this upper 
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Figure 2: The curve on the right specializes Eq. m for Q = 2.5, r+ = 100 , and /3 = 697/3. The curve on the left results from 
dividing the previous one by r, as given by Eq. ®. 



Figure 3: The first two curves from left to right, with A = 8 and 10, respectively, and ro = 25/70, depict the ratio pm!px 
described by Eq. CQ}. Though not apparent in the figure these ratios do not vanish asymptotically, they tend to r_ = 
(25/70) e“^, and r_ = (25/70) e“®, respectively. The third curve, shows, for comparison, the evolution of the ratio for the 
ACDM model, namely, pmjpK = ro a~^ with an identical value for ro. Again, the interaction is seen to alleviate the coincidence 
problem. 


bound on A also ensures Q > 0. 







Inserting m in o and assuming re = —1, we get 
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Q= 3 - 


Aa 


(1 + a)^ 


H 


1 -|- T— 


whence the corresponding differential equation for the DM density is 


«Pm + > 3 + 


A a 


{l + ay 


- 3 


1 


1 -j- 7-_ 


Pm — 0 ■ 


( 11 ) 


( 12 ) 


This has analytical solutions in the limiting cases, a —?> 0 and a oo; namely, pm oc a“3r+/(i+T-+) ^ 

a-3»'-/(i+»'-)^ respectively, where r_ = roe~^^^ > 0; i.e., r never vanishes. Notice that, as in the previous case, the 
limiting solutions are rather reasonable; for small a, only a tiny amount of energy can be transferred to DM whence 
Pm decreased practically as though Q vanished, pm(a <C 1) ~ a~^. For large a, a large amount of energy has already 
been transferred, significantly compensating the “natural” decrease of density; hence, Pm(a 1) — constant. This is 
apparent in Fig. 4, which depicts the numerical solution (IT^ for the cases A = 8 (i.e., r_ = (25/70) e and A = 10 
[i.e., r_ = (25/70) e“®], respectively. 



Figure 4: Left panel: Numerical solution to Eq. m in terms of the dimensionless quantity pm = Kpm/iiHQ) with k = SttG, 
for the choice A = 8. Also shown is the curve of px = pmir. Right panel: The same but for the choice A = 10. 


C. Third parametrization 


Similar to the previous one is 


r = r_ exp 



(13) 


where A > 0 must be lower than ~ 10.125 for Q to be positive and alleviate the coincidence problem. Here, 
r_ = rg Note that r+ cannot be arbitrarily large because it fulfills r+ = rg . 


Figure O depicts the ratio in terms of the scale factor for A = 8 and A = 10 as well as the corresponding plot for the 
ACDM model. 

The interaction term in this case, with w = —1, reads 


Q 



TT __ 

1 + r_ ■ 


(14) 


Thus, the differential equation that governs the evolution of the DM energy density , assuming w 


a Pm 



2 Aa 
(l + a)3 


1 


1 -(- f_ e^/(i+“)^ 


= 0 . 


— 1, reads 


(15) 


In the limit situations a ^ 0 and a —>• oo, one has pm oc a 3’'+/(i+’'+) and Pm oc a 3r'_/(i+r_)^ respectively. A plot of 
the numerical solution for A = 9.89 is shown in Fig. [51 
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Figure 5: Plot of the ratio pm!Px as a function of the scale factor as given by Eq. (1131) . The first two curves from left to right 
are for A = 8, A = 10, respectively, with ro = 25/70. Also shown for comparison is the analogous curve for the ACDM model 
with the same value for ro. 



Figure 6: Numerical solution of Eq. GSl) in terms of the dimensionless quantity = K,pm/{HI q) with k = SttG for the choice 
A = 9.89 and r_ = 0.031334. Also shown is the corresponding curve of px- 


III. CONSTRAINING THE FREE PARAMETERS 

Thus far, we did not resort to any specific theory of gravity. We have just used the FRW metric and the conser¬ 
vation equations (l)-(3). Accordingly, the parameters entering the above parametrizations have been left largely 
unconstrained. In order to most effectively restrict them one should use some gravity theory alongside the full set of 








observational data of supernovae type la, baryon acoustic oscillations, history of the Hubble factor, cosmic microwave 
background, growth factor, and so on. However, this is beyond the scope of our work as we only aim to illustrate how 
reasonable expressions for Q can be obtained from the qualitative evolution of the ratio Pm!Px- This is why we shall, 
for simplicity, assume tc = — 1 throughout and limit ourselves to quickly find acceptable values for the remaining 
parameters. This will be done by using Einstein gravity and the redshift at which the transition from decelerated to 
accelerated expansion took place. 

Recently, Farooq and Ratra [l^, based on 28 independent measurements of H{z), determined with unprecedented 
accuracy that redshift, namely, Zda = 0.74 ± 0.05 (Icr). The corresponding scale factor, ada = (1 + Zda)~^, lies in the 
interval 0.56 < ada < 0.59, being 0.57 the central value. 

From the deceleration parameter, q = —ii/aH^, and the Einstein field equations for a spatially flat universe, 


='^P, and - = (p + 3p), 
d ad 

where p = pt + Pm + Px, P = Px = wpx, we obtain 


9 


1 

2 


[(Pbo/Pmo) + l]r + 1 + 3w 

[(Pho/Pmo) + l]r + 1 


(16) 


(17) 


In the limiting cases r 1 (i.e., a <C 1, matter-dominated epoch) and r ^ 1 (a 1, far future), one follows q = 1/2 
and q ^ {1 + 3w)/2, respectively, as it should. In the particular case w = —1, it reduces to 



1.2r - 2 

F27TT 


(18) 


where we have considered Pbojpmo = 0.2. 

The latter expression of q together with the above value of ada can be employed to get acceptable values for the free 
parameters in Eqs. dSD, (HOD, and (m. 

For the first parametrization Eq. ®, Fig. [3 depicts the dependence imposed on the parameters /3 and a by the 
observational constraint q{ada) = 0. For simplicity we have taken r_ = 0. 



Figure 7: The parameter fi in terms of a as imposed by the constraint q{a = 0.57) = 0 on the first parametrization, Eq. (f^. 
with r_ = 0. 
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Then, the deceleration parameter vanishes at Qda = 0.57 for any point in the graph. Since we are dealing with the 
fixing of two parameters a further fixed point is necessary. We take the product to Hq where to and Hq stand for 
the age of the Universe and the Hubble constant, respectively. Measurements of the age of globular clusters and 
other old objects place the former in the interval 12 < <o ^ 14 Gyr M- As for the latter, it is safe to say after the 
results of the Planck Collaboration [2^, Hq = 67.3 ± 1.2 Km/s/Mpc, and the local measurements of Riess et al [2l|, 
Hq = 73.8 ± 2.4 Km/s/Mpc, that the said product must lie in the range 0.81 ^toHo ^ 1.09. Here, 


to Hq — 


da 


0 a ^ 11^0 ^ ^ U flmO/(^) 4” (1 bo Umo) 


(19) 


where /(a) = pm/po is given by the right-hand side of Eq. ([9|) with r_ = 0 and r+ = ro(l -|- (3), 
g(a) = Pxjpxo = ?'o/(a)/?') H&o = 0.05, and 12^0 = 0.25. Thus, only that pair of values (a,/3) lying on the 
graph of Fig. [7] such that the corresponding product to Hq lies in the above interval can be considered for the 
parametrization. There is a further qualitative restriction. In order that dark energy does not spoil the standard 
picture of primeval nucleosynthesis it must be present only in a comparatively small amount at that epoch, namely, 
ilx < 0.045 [l^. This implies that /3 must be large, no lower than, say, 50. At any rate, given the absence of a second 
strict observational constraint at the background level, there is a certain latitude in choosing the values of both 
parameters, i.e., properly speaking there is not a “best fit” as a handful pair of values (a,/3) fits the data equally 
well. Thus, only after a substantial improvement on the observational estimation of Hq and to, it will be possible 
to determine the best fit. At any rate, for the sake of illustrating the method, we take a = 2.8 and /3 = 111.33 
(consequently, r+ cs 40.12), and this corresponds to to Ho — 0.98. 


Likewise, for the second parametrization setting q{ada) = 0 leads to A = 11.25, thereby r_ cs 0.0013 and r+ ~ 112. 
Similarly, for the third parametrization one obtains A = 9.89 whence r_ ~ 0.030 and r+ ~ 594. As it can be 
numerically checked, the product to Ho yields 0.99 and 0.98 for the second and third parametrizations, respectively. 
Thus, in both cases it falls within the allowed range. 

Figure [5] collects the evolution of the deceleration parameter for the three parametrizations [Eqs. ([5]), (flUl) . and m, 
respectively] for the above parameter values. Though not shown in the figure, qir^) cs 0.46 for the first one and 
cs 0.49 for the other two. Likewise, q(r-) ~ —1 for the two first parametrizations and cs —0.95 for the third one. 



Figure 8: Evolution of the deceleration parameter in terms of the scale factor for the first, second and third parametrizations. 
In drawing the curves, we used w = —1, ro = 25/70, and pbo/pmo ~ 0.2. 


To establish whether the above parametrizations, alongside the parameter values just obtained by imposing q{ada) = 0, 
are acceptable —at least at a good qualitative level — we next determine H{z) for each parametrization and contrast 






it graphically with the history of the Hubble factor. 

From the expressions q = a/aH^ and H = aja, it follows 
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77 = iJoexp|y [1 + g(a;)] dln(l + a;)| , (20) 

where Hq is Hubble’s constant. Expression (1201) provides us with the Hubble history for every parametrization of 
Pm/Px via the corresponding deceleration parameter. 



Figure 9: Hubble history corresponding to the three parametrizations, Eqs. (m, ng, and m- For the curve corresponding to 
the first parametrization Eq. we have assumed a = 2.5 and P = 697/3, for the second parametrization Eq. (HOI) . A = 8, and 
for the third prametrization Eq. (Ha A = 9.89. In drawing the curves we have taken w = —1, ro = 25/70, and pbo/pmo = 0.2. 


Figure[9]depicts the Hubble history (normalized to Hq) in the redshift interval 0 < z < 3 for the three parametrizations: 
(0), nni), and (ITH)) . They begin to differ only after z = 2. Likewise, Figs. ITUl and ITT] show the ratio H{z)/{1 + z) 
for 0 < z < 24, for the same three parametrizations confronted with the 28 observational H{z) data collected in [l8l| . 
In Fig. [Tni we have used Hq = 68Km/(s ■ Mpc) [1^, a value within Icr of the one reported by Ade et al., namely, 
67.3 ± 1.2 in the same units 1^. In Fig. [TTl we used the central value Hq = 73.8Km/(s ■ Mpc), obtained by Riess 
et al. by local measurements [^. The graph for the ACDM model is not shown because it essentially overlaps the 
other three. 


A. The relative weight of Q 


As mentioned in Sec. H, we do not expect the interaction term Q to be larger than the second term in the conservation 
equation for the DM energy density, Eq. (2). Then, assuming w = —1, every parametrization should fulfill 


3pmH 

Q 


3 (1 + r) 

3 + (ar'/r) 


> 1 . 


( 21 ) 


This is satisfied by the three parametrizations —Eqs. (0, (nni) and m- For the first one. Fig. [12] illustrates (using 
the parameter values derived in the previous section) how the relative importance of Q, i.e., (C~^); increases with 
expansion and is bounded from above to be lower than unity (in fact, in this case, C ^ 15 for the whole cosmic 
history). Notice that the graph is fully consistent with the fact that | p'^ \ is much larger for a <g; 1 than for a around 
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Figure 10: Plot of the ratio H{z)/{1 + z) assuming Hq = 68Km/(s • Mpc). 
second and third parametrizations [Eqs. (Q, (IIOI). and (1131) 1. respectively. 
The H{z) data were borrowed from Refs. 2^ - [291 ]. 


Black, blue, and red lines correspond to the first. 
The data error bars indicate la confidence level. 


and beyond unity (see Fig. n. 

A rather similar situation occurs for the second and third parametrizations [Eqs. (fTOl) and dT^ j. This is why we omit 
the corresponding graphs of C- 


IV. TWO “POSSIBLE” DEGENERACIES 

In principle, any given ratio r{a) may be equally associated to two different scenarios; viz., the one considered so far 
(a nonvanishing interaction term Q with constant EOS w for the DE component) and also to the case of vanishing 
Q (no interaction) but with w given by a suitable EOS, different from the previous one. Therefore, one could not 
tell whether such a ratio points to the existence of an interaction or, perhaps, to a more involved EOS. In other 
words, a possible degeneracy between Q and w arises. Here, we shall argue that this degeneracy is only apparent 
because it can be readily broken at the background level since observable quantities, such as the Hubble factor and 
the deceleration parameter, significantly differ from one scenario to the other. 

Indeed, by setting Q = 0 in the set of equations (I)— ([3]), the EOS that produces the same r{a) as the one associated 
to a nonvanishing Q is easily found. 


w{a) 


r' o' 

' rm 

r Pm 


( 22 ) 


For illustrative purposes, let us focus on the second parametrization of r, namely, m, and the corresponding 
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Figure 11: The same as Fig. fTOlbut using Hq = 73.8Km/(s • Mpc). 


interaction term is given by (HU. With the help of (fT^ we obtain 


w{a) 


(l + a)2 


3 + 


A a 


(l+a)^ 


- 3 


1 -|- T— (i+“) 


(23) 


Figure [13] depicts the dependence of w on the scale factor. This variable EOS produces the same r(a) as the case 
with interaction, Eq. (ITpl. However, this EOS seems to be at variance with observation because it contrasts with the 
findings of Serra et al. [^, who using a variety of observational data conclude that the EOS, varies, if at all, very 
little from z = 1 (i.e., a = 0.5) up to the present. 


Let us consider the Hubble factor for the scenario with (5 = 0 and the EOS 

H{a) ' 




= \ flboa 3 3 + ilxo exp 


-3 


1 + w{x) 


dx 


(24) 


where flbo = 0.05, flmo = 0.25, and i2x0 = 0.70 are the current fractional densities of the various energy components. 


Figure [HI depicts the evolution of the Hubble factor associated to Eq. (1^ (solid line), and the one corresponding to 
interacting scenario with w = — 1 (dashed line), both are normalized to Hq. Both graphs depart more and more from 
each other as the scale factor decreases. Further, since ff(a) significantly differs between both models, the luminosity 
distance is also affected, as well as all other cosmological observable quantities (e.g., supernovae la brightness, age of 
the Universe, background acoustic oscillations, etc.). So, the degeneracy gets already broken at the background level. 
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Figure 12: Evolution of the ratio (1211) with expansion for the first parametrization, Eq. dSJ. The larger the scale factor, the 
smaller C. In drawing the curve we have used ro = 25/70, /? = 111.33, and a = 2.8. 



a 


Figure 13: Evolution of the EOS corresponding to the second parametrization, Eq. (O, assuming Q = 0. In drawing the 
curve we have taken A = 11.25 and r_ = 0.0013 (the best fit-values found above; see text). 


The consideration of the deceleration parameter strengthens this point further. Equation dm applies to both 
instances, viz., the model with Q — 0, and with r and w given by Eqs. (IT0|) and (l23)l . respectively, and the model 
with Q given by Eq. dm, where w is fixed to —1. The observable q significantly varies from one model to the other 
because w(a), which enters dm, does also, while all the other quantities in dm are identical. Figure ITS] makes this 
explicit. The difference is already significat, though not large, for a = 0.57 (the value observationally determined by 
Farooq and Ratra [l^ at which q vanishes). For the Q = 0 model, q{a = 0.57) ~ 0.046, while for the other model 
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Figure 14: Evolution of the Hubble factor normalized to its present value, Hq. The solid line corresponds to Eq. (HU with 
w{a) given by (113). The dashed line depicts the Hubble factor associated to the case of the same r(a) but with Q given by HB) 
and w = —1. In drawing both graphs we assumed H(,o = 0.05, HmO = 0.25, = 1 — Hbo — HmO, A = 11.25, and r_ = 0.0013. 



Figure 15: Evolution of the deceleration parameter for the model with <5 = 0 (solid line) and for the model with Q given by 
dill) (dashed line). Both models share the same ratio r, equation dlOl) . but their respective EOS are different -see text. In 
drawing the graphs we took p6o/PmO = 0.2, A = 11.25, and r_ = 0.0013. 


q{a = 0.57) ~ 0.0031, with the former figure being about 15 times the latter. Again, the degeneracy is readily broken 
at the background level. 

Here, we have just considered for illustrative purposes the interaction term associated to the second parametrization. 
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An identical conclusion can be attained for the other interaction terms above. 


Another apparent degeneracy is as follows. In principle, one could think that it would be possible to reproduce the 
same expansion history H{z) and, consequently, the same deceleration parameter q{z) in a noninteracting scenario 
{Q = 0) by using a suitably tailored EOS w{z), for the dark energy component in an interacting scenario (Q ^ 0) 
with w = constant. Obviously, now the ratio riz') will differ from one scenario to the other. Since H{z) is the same 
in both, they will share the same luminosity and angular distances whence the models will be indistinguishable at 
the background level. 

However, as we will see, this possibility does not exist in general, nor does it exist in the particular cases of the three 
parametrizations above or for parametrizations dSO]) and (1321) considered below. Indeed, by equating the right-hand 
sides of Friedmann’s equation for spatially flat universes corresponding to the interacting model, and assuming re = — 1 
(for simplicity), with the nointeracting one Q = 0 but with an undetermined EOS w(a), it follows that the latter 
quantity is given implicitly by the expression 

/ ^ da = i In {ro [/(a) - a”^] -b g(a)} , (25) 

where /(a) stands for the DM ratio Pm{a)/Pmo in the interacting model [e.g. by the right-hand side of Eq. ([9]) in the 
case of the first parametrization] and g{a) = Px/Pxo = ^o/(n)/^(n). In writing Eq. (I25L we bore in mind that the 
current fractional densities of the three components coincide in both models. 


Because of at any a < 1, the amount of DM in the interacting model is necessarily less than in the case of vanishing Q 
(in which instance the DM density would obey pm = Pmo a~^), one follows /(a) < a“^, and the curly brackets on the 
right hand-side of (l25l) can be negative in some interval of the scale factor when a < 1. This severely constrains w{a) 
since the said brackets must be non-negative for any value of a if w{a) is to exist at all. Accordingly, if the function 

X(a) = /(a) - a~^ + (26) 

r(a) 

takes negative values in some interval of the scale factor, there will be no w{a) capable of reproducing the H{a) 
function of a given interacting model. This is the case of the three parametrizations ([5]), (fTOll . and (fT^ . as Figs. fTHl 
and [m illustrate. 

As Fig. [TH] shows, the same is true for the parametrizations given by Eqs. (1501) and (1501) (see below). These 
correspond to the widely used interaction terms Q = ^Hpm and Q = ^Hpx, respectively. Here, ^ is a semipositive, 
small, constant parameter. 


In general, it is reasonable to expect that for any Q > 0 interaction the /(a) function may be “piecewise” expressed as 
/(a) = with 0 < e -C 1 a constant parameter in every small-scale factor interval but such that it slightly varies 

between two neighboring intervals; the smaller the interval, the better the approximation is. Then, the expression 


X(a) ~ a 




(27) 


will be valid in every small interval of the scale factor. Accordingly, bearing in mind that dr(a)/da < 0, and that 
for r > (0.25/0.70)^/^ ~ 0.71, DM dominates over DE (otherwise, the large-scale structure we observe today could 
not be accounted for), it follows that r{a) > 1 for a small enough. Then, for some scale factor, say a* < 1, 
x(a) will be negative in the interval 0 < a < a*, and no EOS w{a) will exist such that the corresponding non¬ 
interacting model can reproduce the same history of the Hubble function of a given interacting model with constant w. 


On the other hand, even wildly assuming the existence of such EOS, both kinds of models will be readily told apart 
at the perturbative level without the need to perform a detailed analysis. Indeed, because in the interacting model 
DM is continuously added to that component from the DE (recall that Q > 0), as mentioned above, the amount of 
dark matter in the past was necessarily less than in noninteracting scenarios at the same redshift. This is illustrated 
by Fig. [T51 where the ratio rj = Pm{z)IPmO (1 + z)^ vs. z is depicted for the case of the second parametrization, Eq. 
in the redshift interval 0 < z < 20. On the other hand, since the present amount of matter (baryonic, ftbo and 
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Figure 16: Evolution of function Xj Eq. (1261) . in terms of the scale factor for the case of the first parametrization, Eq. ((S}. As 
is apparent, x > 0 only for a > 0.4, thereby no EOS of any noninteracting model can reproduce the Hubble history, H{z), of 
the model corresponding to the first parametrization. In plotting the graph we have used ^Ibo ~ 0.05, HmO = 0.25, r_ = 0, 
r+ = 40.12, a = 2.8, and /3 = 111.33. 



a 


Figure 17: Solid line: Evolution of function x, Eq. (1261) . in terms of the scale factor for the case of the second parametrization, 
Eq. (HOI) . Dashed line: The same but for the second parametrization, Eq. H3I) . As is seen, x{o, < 0.32) < 0 and x(a iS 0.30) < 0 
for the second and third parametrizations, respectively. In drawing the graphs we have used the value of the parameter A that 
best fits the observational data [q{a = 0.57) = 0], i.e., 11.25 and 9.89 for the second and third parametrizations, respectively. 
In both cases we have assumed Qbo = 0.05 and = 0.25. 


dark Umo) has to be the same in both models, the growth function, fg = d In Zl-|_/d In a in the interacting scenario 
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a 


Figure 18: Solid line: Evolution of function x, Eq- dsn, in terms of the scale factor for the case of parametrization dMJ- 
Dashed line: The same but for parametrization For both parametrizations x is negative in a considerable interval of the 
scale factor. In drawing the graphs we have considered ro = 0.25/0.70, w = —1, and ^ = 0.01. In both cases we have assumed 
n^o = 0.05 and 12^0 = 0.25. 


must be larger at any redshift than in the ACDM model, while the same function in the noninteracting scenario will 
be, at most, as large as in the latter model. While current measurements of fg{z) are not accurate enough to do the 
job, future measurements —possibly the next generation of data—will, hopefully, be precise enough to discriminate 
between the said scenarios. This may well be achieved by the Euclid mission [^. Thanks to the latter, fg{z) is 
expected to be determined with an accuracy between 1% and 2.5% in the redshift interval 0.5 < z < 2 (see Figs. 
2 and 3 in Di Porto et al. and Figs. 14 and 15 in Amendola et al. [s^). In summary, if, by any chance, an 
interacting model and noninteracting one happen to share an identical expansion history, in order to tell one model 
from the other, it will suffice to see whether the experimentally determined growth function lies above or below the 
theoretical growth function of the ACDM model. If it lies below, the interacting model gets automatically discarded; 
if it lies above it is the noninteracting one that should be disposed of. Altogether, sooner or later observational data 
will constrain fg{z) sufficiently enough so that without explicit calculating the said function for any of the two kind 
of models it will possible to discriminate one from the other with full certainty. 

Likewise, the maximum of the matter power spectrum of interacting models will be shifted to larger scales with 
respect noninteracting models. Indeed, the redshift of matter-radiation equality Zeq determines when subhorizon 
density perturbations start to grow. A model with a larger fraction of DM (i.e., a noninteracting one) will reach this 
equality earlier than a model with a smaller fraction of DM (an interacting model with Q > 0). The power spectrum 
of matter density perturbations has its maximum at fc = 2TT/dH(zeq)i where dniz) is the radius of the horizon at 
redshift z. In any model, when Dm + flh = flrad (1 + Zeq) the perturbations start to grow. Because of the Meszaros 
effect [S^, the lower the fraction of DM, the later the growth of density perturbations will commence, making the 
models differ further. Eventually, data coming from large-scale structure observations will render this distinction 
feasible. 


V. SOME PROPOSED INTERACTION TERMS 

As mentioned in the Introduction, several interaction terms have been proposed in the literature over the years, all 
of them —due to of our lack of knowledge of the dark sector at the fundamental level— were based just on heuristic 
arguments and mathematical simplicity. 
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Figure 19: Evolution of the ratio t] = pm{z)/pmo (1 + z)^, where pm{z) is the dark matter density in the case of the second 
parametrization. The solid line corresponds to A = 10, the dashed one to A = 8. In drawing both curves ro = 25/70 was 
assumed. 


Here we will focus on two sets on interaction terms. The first set is [EHi, 


Q = ^Hpm, Q = iHp^, Q = ^H{pm + Px), Q = ^H , Q =-^(prn + p,y, (28) 

Pm + Px 

where ^, stands for a constant and small, | ^ |<C 1, dimensionless parameter. 

A second set, given by [36 1 

Q ~ ^Pm 5 Q ~ ^Px 5 and Q = r(/?^ + p^'j , (^9) 


where the constant T has the dimensions of inverse of time and was inspired by proposals about the curvaton decay 
and the decay of DM into relativistic particles [s^ . 


By equating the right hand-sides of (|4|) and (l28l ll and integrating the resulting expression, with w = constant, one 
obtains 


(1 -I- c) ro 

1 -I- c -I- ro — ro ’ 


(30) 


where c = 3i(;/^. 


As it can be readily checked, for 0 < ^ <C 1, r is non-negative for any value of the scale factor, see Fig. [501 Further, 
it does not diverge when a —^ 0; in fact r+ = —(1 -f c) > 0. This is apparent in the left panel of Fig.l of Ferreira et 

al. [13. 

For ^ < 0 satisfying | ^ |<C 1, r takes negative values and diverges when a = a, with 


ro 


l€l(i+c) 


Q: 


1 -I- c -I- ro 


< 1 . 


(31) 
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Figure 20: Evolution of pm /for ^ = 0.1 and ^ = 0.01, according to Eq. (I30II . In drawing the curves we have assumed 
u; = — 1 and ro = 25/70. These overlap for a > 0.75. 


Proceeding along parallel lines with the second interaction term, Eq. (I28l 2~). one follows 

_ [(1 + c)ro + 1] - 1 


1 + c 


(32) 


As in the previous case, for 0 < ^ <C 1, r is non-negative for any value of the scale factor, see Fig. By contrast, 



Figure 21: The same as Fig. [20] but for Eq. (1321) . 

for ^ < 0, with I ^ |<C 1, it becomes negative for a > o* = [1 -I- ro(l -f > 1. Further, regardless of the sign 

of r+ diverges (see middle panel of Fig. 1 in 

From Eqs. (jl]) and (E5131. one is led to 


dr 


da 

a 


-I- (2 -I- c)r -I- 1 


( 33 ) 
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Figure 22: Two numerical solutions to (1331) . The lower and upper curves correspond to ^ = 0.1 and ^ = 0.01, respectively. In 
drawing the curves we have taken ro = 25/70 and w = —1. 


Two numerical solutions to this equation are depicted in Fig. 1221 Though not shown, r remains finite in the a —>■ 0 
limit. In fact, = (2/3^) — 1 + \/(974^^y^^'(370, is apparent in the right panel of Fig. 1 of [1^ and in Fig. 2 of 
Olivares et al. [d^. 


For the interaction term given by Eq. (I28l 4l one obtains 


1 + c 


l/c 


= ro ro + 


1 + c 


l/c 


.C(l+c) 


(34) 


Figure [231 shows the behavior of r for ^ = 0.1 and ^ = 0.01. Both graphs essentially overlap; on the other hand, 
there would have been practically no difference if we would have used ^ = —0.1 and ^ = —0.01 instead. In this case 
r+ diverges for any sign of This is consistent with the fact that (128141 reduces to (l28]2l when pm ^ Px', he., when 
a —>■ 0. 


Let us consider now the interaction term (I28l 5l. With the help of the conservation equations (2) and (3), it can be 
cast as 


Q 


S^Hpm 


1 + r + w 
r 


Following the same steps as before, we find 


dr 

^ (1 + r)2 + wr (1 + ^) + 



This equation integrates to 


a = exp [f (r) - F(ro)] , 


(35) 


(36) 


(37) 


where 


F{r) 


2 

3\/| w I \/(l + ?)2 \ w\ -4^ 


tanh ^ 


(1 + 0 I w I -2g(l + r) 

vT^\/(i + 0^ K I -4C 


(38) 
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Figure 23: Evolution of the ratio pmjpx vs the scale factor, according to (1341) . for ^ = 0.1 and ^ = 0.01. In drawing the curves 
we have assumed ro = 25/70 and w = —1. 



Figure 24: Ratio pm/px vs the scale factor according to Eq. (I37II . The upper and lower curves correspond to ^ = 0.1 and 
^ = 0.01. In drawing the curves we have used ro = 25/70 and w = —1. 


Figure [24] shows the evolution of r(a) for ^ = 0.1 and ^ = 0.01. In this case, r+ remains hnite, and for w = —1, it 
reads 


(1-e) + x/(l + C)^ - 

2 ^ 


(39) 


Obviously, ^ cannot be negative. 

Altogether, the interactions terms given by Eqs. (l28lll. (|28l 3l. and (I28l5l fulfill the criteria of Sec. II provided ^ > 0; 
not so the other two since r_|_ diverges. 


As for the interaction term given by Eq. (12^11. we proceed as follows. Let E = yiLoi where 7 denotes a dimensionless 
constant. Next, we equate the right-hand side of the resulting expression for Q with the corresponding one of (jH). 
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After simplifying and introducing the ansatz H = HQr°‘ (piecewise valid for a = constant in the range 0 < a < 1), 
we arrive at 


dr da 

3wr + 7 ri““(l + r) a 


(40) 


It is seen that dr/da will be negative in general only if 7 < 0. However, as the left panel of Fig. 1 in suggests, 
r+ diverges. 

A similar situation occurs for the interaction terms (1^ 21 and (1^ 31. Thus, we infer that the three expressions for Q 
in (ESI) fail as physically acceptable interaction terms in the dark sector. 


VI. CONCLUDING REMARKS 

Ideally speaking, any valid expression for the interaction between the DM and DE components could be obtained 
by writing the equations of motion after varying the corresponding Lagrangian. However, one cannot propose such 
Lagrangian expression without, implicitly or explicitly, presuming a good deal about the microscopic nature of 
these components. Worse than that, as recently noted, because two fluids should not be considered ideal on their 
own if they interact with one another, it is unclear if Lagrangian (even if is just effective) could be ascribed to 
such a situation (^ . This is why we have followed a diverse route when looking for phenomenological, but useful, 
expressions of the said interaction term. 

Any interaction in the dark sector is bound to change the ratio between the energy densities of DM and DE 
with respect to the very particular case of no interaction, as it is directly connected to it [see Eq. O]. Thus, in 
Sec. H, by parametrizing this ratio on simple, sensible grounds, reasonable expressions for Q follow, as Eqs. 

CD, and m illustrate. Here we wish to stress that this method rests solely on the FRW metric and the energy 
conservation equations and not on any specific cosmological model or theory of gravity. (Obviously, other expressions 
for Q differing from the ones proposed here in Sec. H are obtainable using other possible parametrizations of 
Pm/Px-) Likewise, we have restricted ourselves to the case Q > 0 because this helps alleviate the coincidence problem. 

Useful constraints on the free parameters have been readily obtained (Sec. HI) by imposing the vanishing of the 
deceleration parameter at the transition redshift Zda = 0.74 ± 0.05, as recently determined [l^. Admittedly, the 
latter implies the use of one or another theory of gravity —in the present case, general relativity. At any rate, the 
above parametrizations seem to comfortably reproduce the history of the Hubble factor (see Figs. [10] and HH) which 
suggests that these terms should behave reasonably well at the background level because in FRW universes, H{z) is 
directly connected to the luminosity and angular distances. Therefore, they are expected to make good fits to current 
cosmological data from supernovae type la, baryon acoustic oscillations, the shift of the first acoustic peak of the 
cosmic microwave background temperature spectrum, and gas mass fraction in galaxy clusters. Nevertheless, we do 
not claim that these interaction terms should be taken too seriously since they were proposed mainly on illustrative 
purposes and, ultimately, any cosmological model that assumed whatever interaction term must also fit observations 
at the perturbation level, such as the matter power spectrum and the evolution of the matter growth factor; but this 
lies beyond the scope of this paper and will be the subject of a future research. 

We wish to stress that the r(a) expressions based on the criteria laid down in SEc. H do not necessarily are to 
converge to the ACDM ratio r\ oc a~^ at any limit. This is so because ta diverges very rapidly for a —^ 0, which 
contrasts with the condition that r{a) varies very little or remains constant at very early times. Nevertheless, by 
suitably adjusting their free parameters, they can approximate r\ for a ^ 0.1 with great accuracy. 

As Sec. IV illustrates, the possible degeneracies between the interaction term and the EOS of DE are only 
apparent because they can be readily broken at the background and perturbative levels. Specifically, noninter¬ 
acting models are unable to reproduce the same history of the Hubble function of a interacting model (and vice versa). 

Finally, for the sake of completeness, in Sec. V the ratio Pm!Px associated to a handful of interaction terms found in 
the literature has been analyzed. Some of them happen to satisfy the reasonable criteria laid down in Sec. H. 
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